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We characterize the full family of soliton solutions sitting over a background plane wave and ruled 
by the cubic-quintic nonlinear Schrodinger equation in the regime where a quintic focusing term 
represents a saturation of the cubic defocusing nonlinearity. We discuss existence and properties 
of solitons in terms of catastrophe theory and fully characterize bistability and instabilities of the 
dark-antidark pairs, revealing new mechanisms of decay of antidark solitons. 
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I. INTRODUCTION 

Optical spatial solitons are important for their capa- 
bility to beat diffraction and their potential for engineer- 
ing a variety of optical reconfigurable structures includ- 
ing (and not limited to) couplers, deflectors and logic 
gates. In Kerr media where the paraxial propagation is 
described by the nonlinear Schrodinger (NLS) equation, 
only one soliton solution exists once the parameters (i.e., 
nonlinearity and peak intensity or width) are fixed and 
that solution is stable. However, more general nonlinear 
responses can result into bistability of solitons (strictly 
speaking solitary waves) and/or their instability against 
the growth of weak perturbations. 

In this paper we are interested to investigate such fea- 
tures for solitons sitting on a finite background (i.e., dark- 
like) in the context of the Cubic-Quintic NLS (CQNLS) 
with a defocusing cubic and focusing quintic nonlinear 
response. The importance of such model lies in the fact 
that it constitutes the simplest model for a defocusing 
saturable Kerr effect 0, HJ , whose parameters can be ef- 
fectively measured in a relatively simple way by two- wave 
coupling or Z-scan [1, 0|. 

As far as bistability is concerned, the case of dark 
solitons has been investigated with reference to various 
model including the CQNLS @, UHl- In particular Her- 
mann has shown that dark solitons exhibit bistability of 
the second kind, i.e. characterized by solutions possess- 
ing the same full- width-half-maximum (FWHM), albeit 
possessing different amplitudes (and generally invariants 
of motion). This type of bistability was introduced in 
Refs. [l(| El by Gatz and Hermann, to distinguish it 
from the earlier definition [T^ - [T^ | which implies the ex- 
istence of different solutions possessing the same value 
of one invariant of motion (e.g., the power) for different 
values of the internal parameter, typically the nonlinear 



propagation constant (3. The analysis carried out by Her- 
rmann, however, is limited to stationary solitons, while a 
full family of moving dark solitons can exist. The analy- 
sis is further complicated by the fact that the CQNLS in 
the regime considered here is known to possess coexisting 
antidark (or bright on pedestal) solutions @, HH . The full 
family of dark and antidark solitons, once parametrized 
by the velocity, exhibits intriguing features which have 
been overlooked, and which we discuss below. Further- 
more, whether the full family of the dark-antidark mov- 
ing pairs is stable or not, and which are the instability 
mechanisms is still an open problem. Our systematic in- 
vestigation of these problems provides two answers: (i) it 
shows that the criterium demonstrated by Barashenkov 
[l6| for dark solitons provides the correct exhaustive an- 
swer to the stability problem also for antidark solitons; 
(ii) it clarifies that the decay of antidark solitons can 
follow new scenarios in proper regions of the parameter 
space, rather than always blowing-up as conjectured in 
the previous literature, though collapse is in general al- 
lowed even in 1+1 D because of the high power of the 
focusing term. 

Besides being important per se, the knowledge of the 
dynamics of the whole soliton family of the CQNLS is 
also important in view of recent studies which extend 
the investigation of competin g no nlinearities to the non- 
paraxial [l7| and nonlocal fl8l Il9l| regimes. Moreover, the 
full characterization of the soliton solutions and their in- 
stabilities constitute the starting ground for describing 
the feature of dispersive shock waves (DSW, involving 
multiple solitons in the weakly dispersive regime) [20| . 
an active area of research where succesful experiments 
have been recently performed in non-Kerr media under 
different excitation conditions (2l| - [23| . In this respect, 
here we provide the first prediction of a dispersive shock 
wave produced directly by the decay of a solitary wave 
of the CQNLS model. ' 
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II. DARK-ANTIDARK SOLUTIONS 



We start from the (dimensionless) CQNLS equation: 



.du I d 2 
dz 2 dx 2 



which describes a saturable Kerr-like nonlinearity 
through its truncated expansion at second-order in the 
normalized intensity |u| 2 , with a being an external free 
parameter that weights the quintic nonlinear response 
(the smaller a, the weaker the saturation effect). Soli- 
tons of such system have been reported before. How- 
ever, we reformulate the full problem from the begin- 
ning, giving novel analytical formulas which prove con- 
venient for the purpose of our analysis. Solitons cor- 
respond to trans lation ally invariant solutions of the the 
form u(z, x) = \J p(9) exp [i(b(9) + ifiz], where 9 = x — vz 
and /? = g(po) is the nonlinear phase shift experienced by 
a plane- wave background with intensity po = | uo | 2 , in the 
medium where the nonlinear refractive index varies with 
intensity p = \u\ 2 according to the law g(p) = —p+ap 2 /2. 
These solutions depend on two internal parameters which 
we choose, in analogy to general dark soliton solutions 
of the defocusing NLS equation, as po (intensity back- 
ground, which fixes also j3) and v (soliton velocity, which 
fixes also the darkness or brightness of the soliton) . Note 
that, here, the quintic term prevents the simple rescaling 
to po = 1 without rescaling a, and the velocity compli- 
cates further the scenario, so we keep the three param- 
eters free. The modulus p plays the role of equivalent 
" position" , and obeys the standard Hamiltonian dynam- 
ics with "momentum" p = p = dp/d9, 



m . on p 2 



(2) 



V = 2p 



^(p 3 -pl)-(p 2 -p^ + 2k Q (p-p )-^ 
•j Po 



Here Co = vpo and k = ^- + p — ^p\. Once p{9) is 
obtained by solving Eqs. ((2J), the phase profile (f)(9) can 
be found by integrating the following equation: 



(3) 



A. Soliton solutions 

Soliton solutions sitting on the plane- wave background 
p = Po correspond to homoclinic separatrix trajectories 
of Eqs. (J2|), characterized by the energy level T~L = E with 
E = V(po) = — 2cq. Such separatrices emanate from the 
saddle point (p,p) — (po,0) of the Hamiltonian H(p,p). 
For a ^ 0, the potential V(p) — E is a double well cor- 
responding to a double-loop separatrix as shown in Fig. 
[IJa,b). Therefore one has, in general, a coexisting pair of 



dark and antidark solitons that corresponds to the mo- 
tion along the left well p m < p < po (dark solitons), and 
the right well po < p < p a (antidark solitons), respec- 
tively. Here p m < po < p a are the roots of V(p) — E 
(explicit expressions of p m , p a are reported in Appendix 
A). In terms of such roots, we derive (see Appendix A) 
the following explicit solutions for dark solitons: 

ra\ Pm + r Pa tanh 2 [w(9 - O )] 
Pd(. d ) = ——. — ; — ; 2 r Ta — TTTi — ' W 



where r 



l + rtanh 2 [w(6)-l9o)] ' 
(po - Pm)/(Pa - Po) and w 



\J a(p a — po)(po — Pm)l^> is the inverse soliton width, 
while the minimum intensity (dip) is given by the root 
p = Pm- Similarly for anti-dark solitons, we obtain 

la\ P a + rP™ tanh2 i W ( 6 ~ °o)} /,x 

1 + i tanh [w(9 - 9 )\ 
where r and w are the same as for dark solitons. 
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FIG. 1. (Color online) Dark-antidark pair sitting on the unit 
background po = 1, for a = 0.6, v = 0.3: (a) potential V(p) — 
E vs. p; (b) phase-space picture (contour lines of Ti); (c-d) 
Relative intensity (thick black solid line) and phase (thin red 
solid line) profiles of dark (c) and antidark (d) solitons. 

The expressions in Eqs. (jl}-© allow us to obtain, 
by integrating Eq. ([3]) , the nonlinear phase associated to 
the two soliton families. We obtain for dark and antidark 
solitons, respectively 



bd = stan 

w 



b a = — s tan 
w 



' r J^. tanh[ti>(0 - o ] 



'-^tanhK0-0 o )] 
rp a 



0o, (6) 
0o, (7) 



where s 



(Pa-P())(pa-p m ) 



From Eqs. ([B][7|) one can 



easily calculate the phase jump A$ = $(+oo) — oo) 
across the soliton. 

Interestingly, the existence domain of the soliton pairs 
can be described with the aid of catastrophe theory 
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FIG. 2. (Color online) Cusp catastrophe picture for dark- 
antidark soliton pairs. Each curve shows the evolution of 
the parameter a and 6 of the normal form potential V(y) = 
y 4 /4 + ay 2 /2 + by, calculated for a dark-antidark soliton pair 
with fixed internal parameters v and po, and a changing from 
zero up to its critical value a c , where all the curves arrive 
tangentially on the cusp curve (Eq. (|8]), red solid line) that 
bounds the soliton existence domain. 




FIG. 3. (Color online) Potential and phase-space picture for 
v = 0, po = 1, and different values of a: (a-b) a — (ideal 
Kerr case); (c-d) a = 0.5; (e-f) a — 0.9 (close to critical value 
a c = 1, where the right branch of the separatrix becomes 
vanishingly small). 



24] . already applied to characterize dark-antidark soli- 
ton pairs in a different context (gap soliton theory [25$ ). 
In fact, the quartic potential V in Eqs. belongs 
to the A + family, and gives rise to the so called cusp 
catastrophe. According to this picture the potential 
V(p) in Eq. ©, which is of the general form V(p) — 



C4/3 4 + C3/3 3 + C2P 2 + cip, can be cast into the canonical 
|24| V(y) — y 4 /4: + ay 2 /2 + by, by means of the 

- 1/4 y 



form 

change of variable p = (4c i )~ 1/4 y — c 3 /(4c 4 ). Then, in 
the control parameter plane (a, b) (explicit expressions 
of a and b as a function of a, po,v are cumbersome but 
can be easily derived), solitons exist in the inner region 



bounded by the curve (so called bifurcation set [24 

2 

= 0, (8) 



(i) 3 



shown in Fig. [2J Such curve marks the values where the 
critical points (d y V = 0) of the potential become doubly 
degenerate (dyV = 0), and exhibits the characteristic 
shape of a cusp in the origin (three- fold degenerate point, 
dyV = 0). In terms of the original parameter po,v,a, the 
existence condition requires a < a c , with the following 
critical value of the quintic coefficient a c : 



1 

Po 



1 



V 

Po 



(9) 



Taking fixed internal parameters po and v, while chang- 
ing a continuosly, makes the control parameters a and b 
calculated for the soliton to span a smooth curve in the 
control parameter plane (a, b), until at a — a c , the curve 
hits (arriving tangentially) the boundary set by the cusp 
curve [Eq. ©]. Different values of po and v result into 
different control curves, as displayed in Fig. [2] We point 
out that a similar behavior occurs by varying po or v, 
keeping the other two parameters fixed. In particular, 
in the latter case, the existence domain turns out to be 
—v c < v < v c , with the cut-off velocity v c = \fp~a 
obtained by expressing Eq. ^ in terms of v = v c 



UPo, 
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FIG. 4. (Color online) (a) dip intensity (darkness) of dark 
solitons and (b) maximum intensity of antidark solitons as a 
function of velocity v for different values of a. Here the back- 
ground is po — 1. The dashed curves represent the existence 
threshold set by Eq. (0. Note the vertical log scale in (b). 

As an example, we show in Fig. [3] how the typical 
phase plane (potential) changes when a is varied be- 
tween zero and the critical value a c . In this case we 
choose po — 1 and still solitons, viz. v = 0, yielding 
p m = which means that the dark soliton is black re- 
gardless of the value of a, while the antidark is char- 
acterized by a peak intensity p a = 3/a — 2pg. In the 
limit a — shown in Fig. [3ja,b), which represents the 
ideal Kerr case, the potential is cubic, and the separa- 
trix has only one branch corresponding to the well-known 
black soliton solution (p = tanh 2 (x)) of the NLS equa- 
tion, whereas for p > 1 the motion is unbounded and 
no coexisting antidark solutions do exist. In fact the un- 
bounded motion for p > 1 can be thought of as the mo- 
tion in the right well that, however, becomes infinitely 
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deep [V(8 = oo) — > — oo)] and wide (since p a —> oo). 
Viceversa, as a grows from zero, the behavior of the po- 
tential at 9 — oo is inverted, and the right well becomes 
finite, allowing for a eight-shaped separatrix correspond- 
ing to the dark-antidark pair [see Fig. [3jc-d)]. For small 
values of a the antidark soliton has high peak intensity p a 
above the background po, which, however, decreases as 
the saturation parameter a increases. For a approaching 
its critical value a c antidark solitons become shallow [see 
Fig. [3[e-f)], until they reduce to the plane wave exactly 
at a = a c , where p a — > pq. The behavior of dark-antidark 
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FIG. 5. (Color online) Color level plot of FWHM [9 FW hm in 
Eq. (|10|) ] in the parameter plane (a, po) for v = 0.1, dark (a) 
and antidark (b) solitons. (c,d) dark solitons with v = 0.4 (c), 
and v = 0.8 (d). The dashed lines give the critical condition 
a c (p ,v). 

solitons with velocity exhibits intriguing features, which 
can be gathered by plotting the minimum (dip) intensity 
p m of dark solitons (the larger p m , the lower the dark- 
ness) and the peak intensity p a of antidark solitons (the 
higher p a , the brighter the antidark) versus v at constant 
po, for different values of a, as displayed in Fig. |4j For 
non-zero but small velocities, the picture remains quali- 
tatively unchanged with respect to the case v = 0, in the 
sense that the peak intensity p a of antidark solitons de- 
creases continuously from infinity (at a = 0) to p a = po 
at the critical value a c , such that the solitons become 
infinitely shallow (i.e. they reduce to a pure plane wave). 
In this case, however, the dip intensity of dark solitons 
is no longer zero, i.e. they becomes gray solitons with 
darkness po — Pm- The darkness decreases for growing 
velocities v up to a minimum value at the bound veloc- 
ity v c = y/ pa — apQ (obtained by solving Eq. © with 
respect to v for fixed a). This is clearly shown in Fig. 
2J where we summarize the result for a fixed background 
p a = 1. Note from Fig. @£a) that Kerr (NLS) dark soli- 
tons (a = 0) are always darker than the corresponding 



CQNLS solitons of the same velocity (the curve p m {v) 
for a — is always below the other curves relative to 
a^0), and have also a larger phase jump A$ than their 
CQNLS counterparts. 

Interestingly, however, as the velocity grows large 
enough (above v — 0.5 in Fig. it turns out that 
dark solitons can become infinitely shallow at the cut-off 
condition for their existence (p a — > p and hence dark- 
ness tends to zero). Conversely, under the same condi- 
tions, antidark solitons cease to become infinitely shal- 
low, rather reaching a finite minimum peak intensity at 
the cut-off condition for their existence [see Fig. 2Jb)]. 
From Fig. 0Jb) it is also clear that, for small a the bright- 
ness of antidark solitons is nearly independent on the 
velocity. This change of behavior at the cut-off condi- 
tion for the existence (from infinitely shallow antidark to 
infinitely shallow dark solitons) depends the background 
Pq. It can be shown to occur at the value of oiq = 3/(4po)i 
in correspondence of the velocity vq = ^/po/2 (in Fig. |4l 
cto = 0.75 and vq — 0.5). 

The coexistence of dark and antidark solitons consti- 
tute a bistable mechanism such that two different so- 
lutions exist with same parameters (po and v) and dif- 
ferent renormalized invariants M, H, P (see Appendix A 
for their definition). However these soliton families can 
be bistable also according to the definition by Gatz and 
Herrmann, i.e. for fixed a different solutions of the same 
width can exist although they sit on a different back- 
ground pq. This type of bistability was investigated for 
still (v = 0) dark solitons 7}. In order to generalize this 
result to the full family (any v, and antidark case), we 
have calculated the FWHM as 



Ofwhm = — tanh 1 y/ f(p , p m , p a ), (10) 

where f(p ,p m ,p a ) = 2pa P ° po P ^ Pm for dark solitons 
(FWHM taken at half the intensity between p and the 
dip p m ) and f(p , p m ,p a ) = Pa +p~-2p m for antidark soli- 
tons (FWHM at half the intensity between p and the 
peak pa), respectively. The results obtained by mapping 
Ofwhm in the plane (a,po) are summarized in Fig. [S] 
As shown, dark solitons exhibit bistability for sufficiently 
large a, regardless of their velocity. However the range 
of values of a where bistability occurs is greatly reduced 
for large velocities [see Fig. EJc)]. Conversely, as shown 
in Fig. \5ld), antidark solitons are never bistable in the 
same sense (no folding of the level curves of 9pwHM 1S 
ever observed at any velocity). 



III. INSTABILITY SCENARIOS 

Having characterized the features of dark-antidark soli- 
ton pairs, we discuss their stability. We proceed by ap- 
plying a known stability criterium according to which the 
stability is related to the derivative of the invariant mo- 
mentum M (see Appendix A) of the soliton against its 
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FIG. 6. (Color online) (a) Renormalized momentum M vs. 
soliton velocity v, and (b) Hamiltonian H(v) vs. M(v), for 
dark and antidark solitons with fixed po — 1 and a = 0.7. 
The red bullets mark the marginal condition d v M — for 
antidark solitons. 
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FIG. 7. (Color online) Stability maps for antidark solitons 
in the parameter plane (v,a) for (a) po = 1; (a) po = 2. 
The light (cyan) and dark (red) shaded domains correspond 
to unstable and stable solutions, respectively. The dot-dashed 
line sets the value an below which antidark solitons have finite 
brightness at cut-off. The dynamics illustrated in Fig. 8 and 
Figs. 9-10 are relative to the sampled values marked by bullets 
and stars in (a), respectively. 

velocity v. The marginal condition 



separates stable solutions (d v M < 0) from unstable ones 
(d v M > 0). Such criterium, proved by Barashenkov [l6[ 
for dark solitons of the generalized NLS equation is found 
to account also for the instabilities of antidark solitons. 
We emphasize, however, that such criterium accounts for 
real eigenvalues crossing into the right-half plane (usually 
bifurcating from zero eigenvalues associated with neutral 
modes, i.e. symmetries of the model). Other instabil- 
ity mechanisms such as oscillatory instabilities resulting 
from edge bifurcations leading to pairs of complex con- 
jugate eigenvalues entering the right-half plane, need an 
independent characterization. However, we have found 
(numerically) no evidence for instabilities of such kind. 
Moreover, the condition for the existence of solitons [Eq. 
©] turns out to coincide with the region where the back- 
ground plane-wave is modulationally stable. Therefore 
the condition ([TT]) turns out to be an exhaustive criterium 
to assess the linear stability of all the solutions presented 
so far. 

The calculation of the function M(v) for dark soliton 
family shows that such function has always a negative 



slope. Therefore the whole family of dark solitons is sta- 
ble in its existence domain, a conclusion that is fully sup- 
ported by our numerical simulation of the propagation. 

Viceversa instabilities take place for antidark solitons. 
In particular, for small a, it turns out that they are al- 
ways unstable in the whole domain of existence since 
M(v) exhibits always positive slope. Conversely at suf- 
ficiently large a, the momentum M(v) changes its slope 
near the cut-off value for existence v c . This is shown in 
Fig. |51 where we compare the momentum M(v) for dark 
and antidark solitons at a — 0.7 and po = 1. Note that 
the change of slope of the momentum means that two 
antidark solutions with same momentum and different 
velocity exist. These two solutions differ by their Hamil- 
tonian H (see Appendix for its expression). Indeed, as 
shown in Fig. HJb) the Hamiltonian H as a function of 
M is folded, showing an upper (unstable) branch and a 
lower (stable) branch. 
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FIG. 8. (Color online) Dynamics of antidark solitons with 
po = 1, a — 0.7, exhibiting different behavior depending on 
the initial velocity v [see bullets in Fig. 7(a)]: (a) v — 0.5, 
stable propagation; (b) v = 0.4, decay into a nearby stable 
antidark soliton; (c) v = 0, decay into a pair of antidark 
shallow solitons with opposite velocities. 

The velocity which gives the marginal stability {d v M = 
0) clearly depends on po and a. The results of Fig. [5] 
could be repeated for different values of po and a, and 
summarized by drawing stability maps in the plane v,a 
at constant p . These are displayed in Fig. [7] for two 
different values of the background po- As shown a rela- 
tively small island of stability is found in the vicinity of 
the boundary for existence, where they have small bright- 
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FIG. 9. (Color online) Decay of an unstable antidark solitons 
(po = 1, v = 0, a = 0.1) into two symmetric dispersive shock 
fans (trains of dark solitons): (a) color level plot of the in- 
tensity; (b) snapshot of intensity (black solid line) and phase 
(red solid line) at z = 5. For comparison the input intensity 
(renormalized to the maximum of the plot) is shown (solid 
blue line). 
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FIG. 10. (Color online) As in Fig. 7 for v = 0.8. 



ness. In particular stability requires a > a m , where a m 
corresponds to the vertex of the stability island [see Fig. 
[7] . In order to test the validity of the marginal stability 
condition we have made extensive simulations performed 
by means of the well-known split-step method. Such sim- 
ulations confirm indeed that d v M — gives the threshold 
for stability also for antidark solitons. The numerics also 
reveal two basic mechanisms of instability, which are il- 
lustrated below by means of numerical runs performed 
by launching the exact soliton profile, while the pertur- 



bation arises from intrinsic roundoff and discretization 
errors. 

The first scenario is valid for relatively large a, such 
that a stable range of velocities exists. As an exam- 
ple for illustration purpose, we have chosen a = 0.7, 
Po = 1, which yield a range of stable velocities v = 
(0.4644,0.5477) [see Fig. [B[a)]. The dynamics for three 
different values of v, marked by bullets in the map of Fig. 
[7ja), is displayed in Fig. [8j As shown in Fig. [8ja), for 
v = 0.5 the antidark soliton propagates in stable fashion. 
Conversely, when v is decreased just below the threshold 
for instability [v = 0.4 in Fig. Hljb)], the initial soliton 
decays into a different antidark soliton of higher velocity, 
which is stable. However, for lower velocities, such that 
the soliton lies far from the threshold for stability, the 
decay into a pair of stable antidark soliton is observed, 
a process which is accompanied by radiation. Note that 
the instabilities are extremely long-range in both cases, 
i.e. they become manifest only after tens of diffraction 
lengths. As shown in Fig. HJc), the two antidark solitons 
are symmetric in the limit v = in such a way that the 
initial zero momentum is conserved, otherwise they ap- 
pear to be asymmetric (case not shown). The scenario 
illustrated in Fig. [8J holds also for different values of a, 
providing a > a m , i.e. a stable range of velocities ex- 
ists. It is interesting to note that the value po — 3/ (4po), 
which discriminates the fact that antidark solitons can 
be or not infinitely shallow at cut-off, lies in the region 
of stability for any po (see dot-dashed lines in Fig. [7]). 
Therefore antidark solitons with arbitrarily small bright- 
ness that exist for a > o.q around cut-off, are always 
stable, whereas antidark solitons which have finite bright- 
ness at cut-off can propagate stably only for a > a m . 

The scenario discussed above change qualitatively 
when a < a m , where the island of stability shrinks to 
zero. In this case the decay instability towards other an- 
tidark solitons is forbidden because no stable solutions 
exist. In this case, the only stable solutions are dark soli- 
tons and therefore the decay instability occurs towards 
these solutions, even though their shape differ dramati- 
cally (indeed being "opposite") from the input antidark 
shape. Importantly when a is small (weak saturation) 
the antidark solitons possess large amplitude and large 
power P a (or number of particles, see Appendix for its 
definition). Under these conditions we have found that 
the decay instability of the antidark leads to a DSW 
(see Refs. [20l-[23l| and references therein), i.e. an ex- 
panding region filled with fast oscillations which behave 
asympotically as solitons. In the present case, starting 
with a zero-velocity antidark soliton, the decay instabil- 
ity leads to two symmetric DSW fans, as displayed in 
Fig. [S] for a = 0.1. In each of the two fans, the inner 
edge is set by the darkest and slowest soliton, whereas 
on the outer edge the fan is linked to the plane wave 
through a train of solitons with progressively decreasing 
darkness and increasing velocity, which become denser as 
the background is approached. Although this behavior is 
reminiscent of that ruled by the integrable NLS equation 
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(limit a = 0) in the semiclassical regime (i.e., nonlin- 
earity much stronger than diffraction/dispersion) under 
excitation of e.g. a gaussian on pedestal [2lj, it must 
be emphasized that this is, to the best of our knowl- 
edge, the first example where the dispersive shock waves 
occurs directly through the decay of an unstable solitary 
wave of the system. Further characterization of the shock 
fan need to develop the Whitham modulational theory 
for the system, which is beyond the scope of this paper. 
However, simple arguments based on the features of soli- 
tons as shown in Fig. |H let us predict that the fan ruled 
by Eq. (1) is narrower as compared to the one ruled by 
the integrable NLS equation under the same excitation. 
This is due to the fact that the soliton with vanishingly 
small darkness, constituting the outer edge of the fan, 
correspond to progressively reduced velocity as the quin- 
tic nonlinearity grows (i.e., a increases), as clearly shown 
in Fig. Ufa). 

The decay scenario shown in Fig. [§] does not change 
when starting with an antidark soliton with non-zero ve- 
locity (v, M ^ 0), except that the two shock fans become 
asymmetric (the higher the velocity, the higher the asym- 
metry), in both the number of solitons and the velocity 
of the darkest soliton (inner edge of the fan), as shown 
in Fig. [TO! This asymmetry is definitely expected, based 
on the fact that a symmetric configuration would have 
M = and hence would lead to violation of momentum 
conservation. Moreover the asymmetric development of 
the shock is analogous to DSW generated in the inte- 
grable NLS limit when starting from a gray beam with 
non-zero velocity and momentum (23j . 

Finally we point out that, for larger a, yet with a < a c , 
the scenario remains qualitatively unchanged though the 
number of dark solitons generated in the decay of the 
antidark soliton decreases. 

Though the aim of this paper was the full charac- 
terization of solitary waves sitting on a finite back- 
ground, before concluding, we point out that such soli- 
tons can also coexist with bright solitary waves with 
zero pedesta l. Se eking for bright solitons of the form 
u(x, z) = y p(x) exp(if3z), a simple calculation shows 
that the peak intensity of these solutions turns out to 
be p a = 3(1 + \/l + 8f3a/3)/(2a). The existence domain 
of such solutions include arbitrarily small a, and hence 
they coexist with dark-antidark pairs. Given the defo- 
cusing nature of the leading order (cubic) nonlinearity 
this might appear surprising. However, it is not difficult 
to understand that such bright solitons are sustained en- 
tirely by the quintic focusing nonlinearity since p a di- 
verges in the limit a — > 0, where the quintic nonlinearity 
vanishes indeed. Without deepening the study of the 
bright case, already investigated in the framework of the 
CQNLS equation, we expect them to be unstable at least 
in the limit of small a, where they do not affect the decay 
dynamics of antidark solitons or the stable dynamics of 
dark solitons discussed above. 



IV. CONCLUSIONS 



In summary we have discussed the main properties of 
solitary solutions with finite background of the CQNLS 
equation with focusing quintic term. Bistability and in- 
stabilities have been studied for the full family of solu- 
tions obtained in new analytical form, and parametrized 
by the intensity of the background wave and the velocity. 
We have found that the solutions exhibit a non-trivial 
behavior against the velocity such that, depending on 
the value of the quintic nonlinearity, either dark or an- 
tidark solitons, become infinitely shallow at their bound 
for existence. Furthermore dark solitons are bistable for 
any velocity, whereas antidark are never bistable. Finally 
dark solitons are stable against weak perturbations while 
antidark are mostly unstable, exhibiting different mecha- 
nisms of instability. In particular, a novel mechanism in- 
volving the decay of an antidark soliton into a dispersive 
shock wave has been characterized. Further work will be 
devoted to assess how the quintic nonlinearity affects the 
formation and dynamics of dispersive shock waves which 
develop from more general (non-solitary) inputs. Poten- 
tially also the collapse, i.e. blow-up at a finite distance 
(which is known to occur even in the 1+1D case that we 
dealt with, owing to the high order of the focusing term 
[1EH3) cou ld play a substantial role that needs further 
investigation. 

From the experimental point of view, the most natural 
setting for testing these results is the study of parax- 
ial beam evolution in nonlinear optics of centrosymmet- 
ric media, where the quintic term accounts for the sat- 
uration of the Kerr nonlinear index ri2 (ni < 0, defo- 
cusing media), usually quantified in terms of the high- 
order nonlinear index 714 Q, the overall index change 
being An — —\n2\I + n^I 2 , where / is the optical in- 
tensity. While the nonlinear indexes n2,«4 (or equiv- 
alently the nonlinear susceptibilities) depends solely on 
material properties and can be accurately characterized 
by means of consolidated techniques [H, [5[ , the normal- 
ized coefficient a used throughout the paper turns out to 
depend on the input intensity as well, and hence the im- 
pact of the quintic nonlinearity can be tuned by changing 
the optical power [T3|. An other area where the predic- 
tions based on the present model can be relevant and 
can lead to experimental test, is the dynamics of ultra- 
cold atoms (Bose-Einstein condensates) , where the quin- 
tic term arises from higher-order (three-body) atom in- 
teractions [1^, and tuning of the nonlinearities can be 
achieved by means of Feshback resonances. 
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APPENDIX A 



Analytical solutions can be worked out from the equa- 
tion which follows directly from Eq. @ 



/4a 



P = \l 



(12) 



where Q(p) — tz[E — V(p)]. The polynomial Q(p) can 
be expressed in terms of its ordered roots p m , po (double 
root), p a , with ordering p m < p < p a , as Q(p) = (p - 
Pa)(p — Pm)(p — Po) 2 - From Eq. (fl~2"f one obtains the 
following quadrature integral, 



p(8) 



dp 



p(«o) y/Q{p) V 3 



(13) 



where p(0o) = p m or p(9o) = p a for dark or antidark 
solitons, respectively. The integral (|13p gives, upon in- 
version, the solutions given in the text. Explicitly, the 
extremal roots p m ,p a are expressed by: 



3 - 2ap Q - y/{Z - 2ap ) 2 - I2v 2 a 



2a 



= 3 - 2ap + y/(3 - 2ap ) 2 - 12v 2 a 
Pa 2a ' 

Note that, in the limit v = 0, the smaller root p m van- 
ishes, and as a consequence the dark soliton becomes a 
black soliton. 



Once the soliton solutions are known one can easily cal- 
culate the renormalized invariants (momentum, Hamilto- 
nian) which are defined as follows 



M 



! (i' r u ■ //., u ! ) | 1 - — ) dr. ( I") i 



H 



p(x) 



[g(po) - g(p)]dp 



dx, (16) 



and the power or number of particles, P a (antidark) and 
Pd (dark): 



/+00 p-roo 
\u\ 2 ~p dx; P d = po- \u\ 2 dx {17) 

-oo J — oo 



The quantities M and H are those employed in the text 
to assess the stability of soliton solutions. 
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